Fermions in Optical Lattices across Feshbach Resonance 



in 
o 
o 



Roberto B. Diener and Tin-Lun Ho 
Department of Physics, The Ohio State University, Columbus, OH 
(Dated: February 2, 2008) 

We point out that the recent experiments at ETH Q on fermions in optical lattices, where a band 
insulator evolves continuously into states occupying many bands as the system is swept adiabatically 
across Feshbach resonance, have implications on a wide range of fundamental issues in condensed 
matter. We derive the effective Hamiltonian of these systems, obtain expressions for their energies 
and band populations, and point out the increasing quantum entanglement of the ground state 
during the adiabatic sweep. Our results also explains why only specific regions in fc-space can be 
populated after the sweep as found in ref. 1]. 



X3 
O 

o 



> 
m 

in ; 

(N . 

o 
in 

o : 

a ■ 

a : 

o - 
o ■ 



X 
S3 



Ever since the experimental observation of a 
superfluid-insulator transition of bosons in an optical lat- 
tice 0] , there has been an increasing interest in quantum 
gases in optical lattices among different communities. It 
is generally believed that all novel phases in solid state 
as well as many new ones specific to atomic gases can 
be realized in the optical lattice setting. In addition, 
the prospect of using these systems to process quantum 
information has been raised. An exciting new direction 
is the combination of lattice physics with Feshbach res- 
onances, which allow one to turn particles from weakly 
to strongly interacting. Not only will this increase the 
chance of finding new quantum phases, but also increase 
the speed for quantum information processing. 

The first experiment with lattice fermions across a 
Fesbhach resonance has recently been performed by 
Esslinger's group at ETH 0. The authors prepare a 
band insulator in a deep optical lattice with two atoms 
per site (in different spin states) at a magnetic field where 
the scattering length is small and positive. They then 
vary the magnetic field to drive the system adiabatically 
across resonance and find that several bands are partially 
populated near and beyond the resonance. These re- 
sults, simple as they may be, pose many fundamental 
questions. Since the ground state near resonance can- 
not be a band insulator because many bands are par- 
tially populated, is it an insulator or a conductor? If it 
were the former, the band insulator would then be con- 
nected continuously to a new type of insulator yet to be 
identified, and would be only one point in a more intri- 
cate continuous family of insulators. If it were the latter, 
then there must be a quantum phase transition as one 
approaches resonance, which remains to be discovered. 
It is also not clear whether this conductor is a canonical 
Fermi liquid, due to pairing correlations near resonance. 

The experiments at ETH are performed in the weak 
tunneling limit for the lowest band with at most two 
fermions per well on the average. The very weak cou- 
pling limit is of great interest because of its relevance to 
solid state systems and quantum information processing. 
The dynamics in this limit is dominated by the energy 
spectrum of two fermions in a single well. These spec- 
tra have been studied for spherical and cylindrical 



harmonic traps by applying a pseudo-potential method 
to model the interaction, and for spherical traps using 
a two-channel model 0. In this last reference, the pa- 
rameters chosen are confined to the narrow resonance 
regime, while the ETH experimentsP] work with wide 
resonances 0. The studies of ref. in effect deal with 
wide resonances. These solutions, however, are written 
in the relative and center of mass coordinates of the pair 
l^f)^ = |$)r'|x)r' in a well located at R, and are 
inconvenient for the introduction of tunneling between 
wells and for comparing with experiments. In recent 
months, a number of models have also been proposed 
for fermions in optical lattices near a Feshbach resonance 
in the weak tunneling limit. These models describe the 
resonance by a term 6^(R)amt(R)amj.(R-) + h-c which 
converts a pair of fermions in the open channel (the a's) 
in the same (m-th) band at site R into a tightly bound 
closed channel molecule 6(R) (a boson). Any single band 
model, however, automatically excludes many processes 
allowed by symmetry which turn out to be crucial for 
multi-band population near resonance, and fail to de- 
scribe Feshbach resonances in deep lattices at the outset. 

The goal of this paper is to point out the proper Hamil- 
tonian for fermions in deep optical lattices across Fes- 
hbach resonances. Focusing on single well physics, we 
shall seek a formulation convenient for a lattice setting 
and derive exact results for many key properties previ- 
ously unexplored. These results are directly related to 
experiments 1] , force us to rethink many fundamental is- 
sues in quantum many-body physics, bringing out many 
intriguing possibilities. Our formulation also applies to 
bosons and boson-fermion mixtures. 

The Hamiltonian: For deep optical lattices with 
at most two fermions per site, the Hamiltonian is 
H = h-^ -\- T, where /ir is the Hamiltonian of 
two fermions in a deep harmonic well located at site 
R and T is the hopping between wells. Suppressing 
the site index R, the single well Hamiltonian is h = 



mj nj. 



h.c. 



where Em = + * — ^iV^z are the en- 

ergy levels of a harmonic oscillator with frequencies 

{uLlx,^y,^z)^ m = {mx,m.y,mz); alna creates a fermion 
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in the open channel with spin a and energy E^^- The 
corresponding wavefunction is u^{v) = W^^^ y ^ Umiiri), 

Urn (s) 



1/4 -s^/ 



-H„i [t^, where H„-^{x) is the 



Hermite polynomial, and di = ^/fi/Muji. b is the closed 
channel boson with its wavefunction fixed in the ground 
state of the harmonic well 0) ■ is the (unrenormal- 
ized) energy difference between b and the open channel 
fermions, and am,n is the coupling converting a pair of 
open channel fermions into the close channel boson b. 

Since the potential of the optical lattice is V{r) = 
^ y 2 cos(gri), the single particle states on a 

lattice are products of ID Wannier functions in x, 
y, and z. The tunneling term therefore takes the 
form f - -ER,R'E.„..im(R,R')aL.(R)am.(R') + 
ib(R,R')&'''(R-)^(R-')]> where tm(R,R') are the single 
particle hopping energies for band m and tb(R, R') is 
the hopping energy for the closed channel molecule Q. 
In the following, we focus on very weak tunneling (t — )■ 0), 
which already has very rich physics. The effect of tun- 
neling will be discussed elsewhere. 

It is important to understand the functional form of 
<*mn- Recall that in real space the conversion term 
is proportional to / D'I'(R, r)V'|(R + r/2)'0|(R - r/2), 
where D is the tightly bound closed channel pair, 
is the open channel fermion with spin a, R and r 
are the center of mass and relative coordinates of the 
fermion pair. Since D is in the ground state of the 
harmonic well, and since it is of atomic size, we can 
write £)'''(R, r) = J| - Wg ^(-\/2i?i)^(r)&. Decomposing 
both in the harmonic states {um}7 we get Ofm.n = 
23/2^9/4^ n, / As {^/2s)u^l (s), or 

(-l)("'.-".)/^ ^^ m,+n,- + l , 
. n ^d,m,ln,\ 2 ^' 

j—x.y,z V J J J 

when rrij — nj is even, and a is a constant; am,n = if 
any of the {rrij — rij) is odd. Restricting to a single band 
amounts to keeping only the diagonal terms in eq. . It 
will not describe the conversion process properly because 
many processes have been suppressed. 

Renormalization of short range divergences: For 
all resonance models, v contains an infinite constant 
since the actual closed channel pair is approximated by 
a point like boson b. Renormalizing 17 to a finite v* is 
a crucial step for any calculation; the finite parameters 
(z/*,q;) then have to be converted to physical scatter- 
ing quantities such as the scattering length and the 
effective range in free space in order to make con- 
tact with experiments. To illustrate the renormalization 
procedure, we consider cylindrical traps with frequency 
uj.j; = ujy = ti)_L, and define 7 = lOz/uj^, dj_ = yJh/Muj±. 

To perform the renormalization, we start with the 
eigenstate |M/) = + Em,n ^m,„aLT«ii) l^)- The 



Schrodinger equation -ffl^*) = E]"^) then gives 



?7m,n = P 



E~E„ 



E-E„ 



{E - E^^nf 



(2) 
(3) 

(4) 



where i?m,n — E^ + E^. Substituting ^ into Q one 
sees that the sum in eq. (@J is divergent. The divergence 
can be extracted Jjy putting a cut off K* in the sum as 
explained in ref. 9], and we have 



A" 

E-V=-CY,PiK), 

K=0 



F(A-) = 



(5) 



where C 



Mo 



and X 



1 — ^ is the 



excitation energy in the center of mass frame. For large 
K\ the sum in © is EkIo ^i^) = ^ 



where W 



C(l/2) = -1.46. We then have 



K=0 



1 
A'+l 



2VW^ Uil, 

+ %/7C(i), 



E-v* = -CW, ly* = V-2Cy/rK* 



(6) 



With short range divergences removed, eq. ® can be used 
to find the energy levels for given v* and a. 

Relating (i^*,a) to {as, To)- For homogeneous sys- 
tems, the energy of a bound pair —\E\ just below reso- 
nance is given by —aj^ — roK^/2 + k = 0, or \J —E = 

^"TF ~ (^^^) ^' ^^^^'^ l-^l = li^K^/M. This case 
can be recovered from eq. (jSJ) by taking uj±_ with the 
bound state energy held fixed, which means x — > —00. 
Eq.® then reduces to y/—E = A{E — v*), where A = 

2^1^^'2 ■ Comparing it with the free space expres- 
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sion mentioned about, we have the relation 



1 

as 



V2iy*h^ 



y*\ro\M 
2h^ 



(7) 



(notice that the effective range in a two-channel model is 
negative). We can then rewrite eq.JHJ in terms of physical 
parameters and Tq, 



as a± 



(8) 



which enables one to find the energy spectrum for any 
scattering length and effective range. 

Wide resonances: For ""^K, \ro \ can be obtained from 
the method in ref. 6] and is found to be 2.4 nm j^, while 
it is even smaller for ^Li. For typical traps, we then have 
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FIG. 1: Energy spectrum {x — — 1 — 7/2) of two fermions 
in a harmonic well versus inverse scattering length. We have 
chosen parameters used in ref.Q], 7 = 5/6, (see text). 

ro/d± <C 1. This corresponds to the so-called "wide res- 
onances" where the wavelength k^^ of the typical scat- 
tering process {k~^ d± in our case) is much larger than 
the effective range ro- The main characteristic of wide 
resonances is that the system will have very small closed 
channel component /? <C 1 and many properties become 
universal. The former can be seen from eq. I0J, which 
can be summed to [l^ . 

p-^^l + Jl^±Hl,2,,p-.)^-^^ (9) 
V 2 7r|ro| ^ pi 

where $(z,s,x) = y^°° n / ^ is the ^-function. The 
emergence of universality can be seen from (jS)) , which re- 
duces to ^y/2j = — W when \ro\/d±^ <C 1. At resonance, 
Qs = ±00, the energy spectrum is found from W{x) = 0, 
completely independent of any microscopic details of the 
system. Note also that in this limit, eq.® reduces to the 
results of refs. ^ in the corresponding regimes. 

Population of different bands: The probability of 
a fcrmion (with cither spin) in the m-th band is 

PM = 2i.,„,„i» = ^^?!^, (10, 

n n ^ ' ' 

which converges for all values of the energy E. In most 
cases, the calculation can be simplified. For instance, 
for a band with indeces (0,0, m), we get F(0, 0,™) ~ 

where the sum is restricted to values of n such that 
TO -|- n is even Il3^ . 

Our results: (I) Energy spectrum: In figure ^ 
we have plotted the scaled energy x as a function of 
—d±/as using eq.JHJ for parameters of the ETH ex- 
periment: 7 — 5/6 and d± = y^h/ {Mlu±) — 63 nm, 
\ro\/d± — 0.04. Even if we set \ro\/d± — (very 
wide resonance), the results change only by a small 



FIG. 2: Population of energy bands for the process \a) 
\d} \b) in Figure The black solid, dotted, and dashed 
curves correspond to the (0, 0, 0), (0, 0, 1) and (0, 0, 2) bands, 
resp. The grey solid and dashed-dotted lines are for (0, 1,0) 
and (0,2,0). In \a), the (0, 0, 0)-band is fully filled, which 
occupies the entire first Brillouin Zone (black area in in- 
set). In state |6), the bands (0,0,0), (0,0,1), (0,0,2) are 
25%, 50%, 25% occupied. The latter two bands corresponds 
to the grey and dotted areas in the inset. Higher Brillouin 
Zones in a;j/-directions are empty in state |6). 

percentage. The eigenstates (|a),|6),|c) in fig.l) at 
—d±/as = ±oo correspond to the lowest three eigen- 
states of non-interacting fermions affected by the inter- 
action; where \a) = \0)rei\0)cm, \b) = |0, 0, 2)rei|0)cm, 
and|c) = 2-i/2(|2,o,0)^ei + |0,2,0)^ei)|0)cm. States near 
resonance such as \d) in fig.l are more complex combi- 
nations of harmonic states in the relative coordinates. 
Note that \b) and |c) are stretched along z and the xy- 
plane with respect to the ground state \a) in the rela- 
tive coordinate, respectively. It is also easy to see that 
in terms of the single particle basis, both \a) and |6) 
can then be represented as J2m, ^m^,n^ nz), with 
= n^ = (0, 0), where m_L = {mx,my). In particular, 
we have \a) ~ |0, 0), and 

|6) = i|0,2) + yi|l,l) + i|2,0), (11) 

which is an entangled state. Thus, as the ground state 
in ref. 0] evolves adiabatically from \a) — > \d) — > it 
acquires orbital deformation and entanglement. Since the 
difference between \a) and |6) is their wavefunctions along 
z, there will be no change in their momentum distribution 
in kx and ky directions. Their differences will show up 
as a spreading of the original fully filled (0, 0, 0)-band (in 
\a)) to the (0,0, 1) and (0, 0, 2) bands (in state \b)) along 
z as shown in figure 2 [lj| . 

(II) Band populations: In figure 2, we have plotted 
the population of some of the lowest energy bands for the 
branch going from a to 6 at different —d±/as using eq. 
(|10ll . All populations can be worked out analytically in 
the manner explained in [jl. ll2l|.As one approaches reso- 
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nance starting from the band insulator |a), many bands 
are populated. The structure of the entangled state \b) 
is manifested in the 2 to 1 ratio of the population of the 
second and fourth bands m = (0, 0, 1) and (0, 0, 2). Near 
resonance, five bands are significantly populated, which 
make up 85% of the probability. The population of the 
close channel molecule is less then 1% according to 
eq.© for the process \a) — > \b). 

Comparison with the ETH experiments Q: Our 
results show agreement but also discrepancies with the 
findings in ref.[l|. First of all, our results predict that 
the momentum distribution of |6) will differ from that 
of |o) only along k^, with substantial population in the 
(0,0,1) band. These are observed in jlj. However, the 
fraction of particles in the (0,0,l)-band in \b) is about 
20% in ref. 1], where as the exact single well result pre- 
dicts P(0,0,1) = 50%, (see fig. 2). In addition, ref.Q 
did not observe a significant population of the (0, 0, 2) 
band in \b), whereas our results shows P(0, 0, 1) = 25%. 
Recently, we learned from Michael Kohl that the fraction 
(/) of singly occupied sites in ref. 0] may be as large as 
50%. The population of the (0, 0, l)-band according to 
our result will then be (1 - /)P(0, 0, 1) « 25%, which 
agrees with Q within experimental uncertainty of num- 
ber of doubly occupied sites. The absence of significant 
population in the (0, 0, 2) band may be due to large tun- 
neling at that energy, which can lead to particle loss when 
the band is populated. Another discrepancy is that the 
original resonance appears to be shifted in ref. 0| whereas 
it is essentially unshiftcd in the exact single well solution. 
We suspect this "shift" may be due to the imaging pro- 
cess, which first lowers the barrier between wells slowly 
(in order to turn quasi momenta into real momenta) be- 
fore a rapid turn off of the trap. In the unitarity regime, 
where scattering is strongest, there will certainly be a 
re-distribution of particles from higher bands into lower 
ones as the barrier is lowered, which may appear to be a 
shift of the resonance. More experiments, however, are 
needed to clarify the situation. We hope that our results 
will serve as a guide for future investigations. 

Final Remarks: The ETH experimental is yet an- 
other example of the rich and subtle physics of Feshbach 
resonance. While the phenomena may appear benign at 
casual inspection, as explained in the introduction, they 
have profound implications on many-body physics. In 
fact, we are in a lucky situation. Far from resonance, 
fermions in deep lattices can simulate nearly all impor- 
tant models in solid state physics by varying the barrier 
height, and hence the ratio U/t (interaction/tunneling). 
Yet, near resonance, one can have a whole host of dif- 
ferent states whose properties are yet to be determined. 
Resonance physics also makes lattice fermions a fertile 
ground for searching non- Fermi liquids as well as for new 
kind of insulators. Although we have discussed mainly 
about the branch labelled a-b in fig.l, the branch e-a is 
equally fascinating for one can turn a band insulator into 



a Bose superfluid. The prospects are truly exciting. 
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